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The problem of measuring nontrivial static correlations in deeply supercooled liquids made re- 
cently some progress thanks to the introduction of amorphous boundary conditions, in which a set of 
free particles is subject to the effect of a different set of particles frozen into their (low temperature) 
equilibrium positions. In this way, one can study the crossover from nonergodic to ergodic phase, as 
the size of the free region grows and the effect of the confinement fades. Such crossover defines the 
so-called point-to-set correlation length, which has been measured in a spherical geometry, or cavity. 
Here, we make further progress in the study of correlations under amorphous boundary conditions 
by analyzing the equilibrium properties of a glass-forming liquid, confined in a planar ("sandwich") 
geometry. The mobile particles are subject to amorphous boundary conditions with the particles in 
the surrounding walls frozen into their low temperature equilibrium configurations. Compared to 
the cavity, the sandwich geometry has three main advantages: i) the width of the sandwich is decou- 
pled from its longitudinal size, making the thermodynamic limit possible; ii) for very large width, 
the behaviour off a single wall can be studied; iii) we can use "anti-parallel" boundary conditions to 
force a domain wall and measure its excess energy. Our results confirm that amorphous boundary 
conditions are indeed a very useful new tool in the study of static properties of glass-forming liquids, 
but also raise some warning about the fact that not all correlation functions that can be calculated 
in this framework give the same qualitative results. 



I. INTRODUCTION 

The sharp slowdown observed in supercooled liquids at 
low temperatures has long been conceptually connected 
to the buildup of structural (static) correlations. Yet, 
due to the amorphous nature of the excitations, it has 
proved rather difficult to identify them and to measure 
their size. For this reason, dynamical correlations [H4l2j 
were detected much before static ones, and only recently 
were structural correlations unveiled, using novel tech- 
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niques [H-[T3|. 

Among these techniques, numerical simulations with 
Amorphous Boundary Conditionss (ABCs), and the re- 
lated point-to-set correlation length £, have proved very 
fruitful [l4|, [l8l - [23| . Implementing ABCs is simple, at 
least in numeric simulations. Consider a set of mobile and 
another one of frozen particles and let the mobile parti- 
cles evolve under the influence of the frozen ones, eventu- 
ally reaching thermodynamic equilibrium. The simplest 
case is when the frozen particles belong to a single equi- 
librium configuration surrounding a spherical cavity of 
mobile particles, of radius R. It is possible then to de- 
fine an overlap q(R) and to measure the similarity at the 
centre of the sphere between two configurations. The de- 
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can check whether or not the finite-size crossovers of the 
correlation functions turn into bona fide transitions. 

Second, when the two walls are very far from each other 
id — > oo), we can study the decay of the overlap off a 
single wall, as a function of the distance z from the wall. 
This is not strictly impossible in spherical geometry, but 
in that case one could be exposed to spurious curvature 
effects that are absent in the sandwich. 

Third, in the planar geometry we can use different 
amorphous boundary conditions on the two sides of the 
sandwich (Fig. [1] bottom), which is also impossible in 
the cavity. This sort of "anti-parallel" boundary con- 
ditions can be used to force a domain wall in the sys- 
tem, and therefore to measure its excess energy and the 
stiffness exponent 9. These quantities are crucial in any 
phenomenological description of the glass transition, so 
that any new tool able to provide information on these 
quantities may be helpful. 



II. MODEL AND SIMULATION DETAILS 



FIG. 1. Cartoon of the sandwich geometry. In the "paralell" 
(or aa) setup (top) both frozen walls are taken from the same 
equilibrium configuration, while in the "anti-parallel" (a/3) 
case, they come from different configurations (bottom). 



pendence of the overlap on the linear size R of the cavity 
yields the correlation lenght £, defined by the crossover 
at R ~ £ among the values q(R) ~ 1 (almost identical 
configuration) and q(R) ~ (statistically independent 
configurations) . 

The original ABC spherical realization can be gener- 
alized to different geometries, where the frozen particles 
do not necessarily form a closed cavity [2(| HH, Hit ■ I n 
this work we study the case of a planar (or "sandwich" ) 
geometry (see Fig. [TJ. As with the spherical geometry, 
in the sandwich we can calculate a point-to-set length 
by studying the sandwich width beyond which the in- 
ternal mobile particles reach ergodicity. In this respect 
our study aims to verify the results obtained in the cavity 
and test their robustness. In particular, we are interested 
to check whether or not the anomalous nonexponential 
behaviour of the point-to-set correlation function at low 
T observed in the spherical geometry [l9| is also found in 
the sandwich. To check how general is this nonexponen- 
tial behaviour is important becoause it is one of the very 
few qualitative thermodynamic landmarks of the deeply 
supercooled phase. 

But the sandwich also allowd us to study cases that 
are out of reach in the cavity. First, in the sandwich 
the width d and the longitudinal size L are independent 
parameters, so that we can perform the limit limit L — > 
oo while the confinement length keeping d finite. This 
thermodynamic limit is clearly impossible in the cavity. 
This limit is important, as by increasing the number of 
mobile particles at constant degree of confinement, we 



We study the soft-sphere binary mixture [25|], a sim- 
ple model of supercooled liquids widely studied before, 
and in which the point-to-set correlation has been com- 
puted using a cavitv. We use the accelerated Swap Monte 
Carlo algorithm [26[ to thermalize the system at tem- 
peratures as low as possible. We run simulations at 
T = 0.482,0.350,0.246,0.202. The first two temper- 
atures correspond to the high-temperature liquid, the 
third is near the "onset" or "landscape-influenced" tem- 
perature [27l | and the lowest temperature lies in the su- 
percooled regime, in which the landscape is dominated by 
minima of the potential energy rather than saddle points. 

The confined system is generated from configurations 
taken from equilibrated pcriodic-boundary-conditions 
runs. These runs were done with density p = 1 and 
box sizes L = 16 and L = 25.3. At each temperature 
we then chose several (from 16 to 24) configurations and 
artificially froze in their equilibrium positions all but M 
particles contained within a region of the simulation box 
in the shape of a box of size 2d x L 2 (we measure d along 
the z axis). 

In order to keep the density fixed within the region of 
mobile particles, it is a standard practice to place virtual 
walls at the border of such mobile regions. What we do is 
the following: taking configurations of the liquid system, 
we place a hard wall potential enclosing the free parti- 
cles. This destroys translational invariance along the z 
axis, but not along the xy planes, creating a sandwich 
of mobile particles surrounded by two infinite walls of 
frozen liquid. 

The main observable we consider is the infinite time 
limit of the local density-density correlations. More pre- 
cisely, we define the overlap q(z;d), as follows: we par- 
tition the simulation box in many small cubic boxes of 
side t, such that the probability of finding more than one 
particle in a single box is negligible. If rii is the number 
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of particles in box i, then, 

q(z; d) = hm ^(t^oH^o + i)), (I) 

where the sum runs over all boxes that lie on a plane 
parallel to the xy plane at the given distance z from one 
reference wall, AT* is the number of boxes in each of those 
planes, and (...) indicates a thermal average. Normal- 
ization is such that the overlap of two identical config- 
urations is 1 on average, while for totally uncorrelatcd 
configurations q=qo = i 3 — 0.062876. 




III. DIFFERENT STATIC LENGTHSCALES 

In this section we study the the overlap, Eq. [TJ in the 
sandwich geometry described above, in which mobile par- 
ticles are confined within a volume 2dL 2 by two walls 
made of frozen particles (top scheme in Fig. [1} . In our 
description, d is the half-width of the sandwich; we be- 
lieve this is the correct variable to compare our results 
(especially lengthscalcs) with the spherical cavity case. It 
is important to note that both walls are made from par- 
ticles taken from the same equilibrated configuration. 

In general, the overlap is a measure of the nonergodic- 
ity of the mobile part of the sandwich due to the frozen 
boundary conditions. When the overlap is nonzero (more 
precisely: larger than its ergodic value go) it means that 
the phase space available to the particles' relaxation is 
reduced by the confinement. It is therefore natural to 
ask how "far away" the walls need to be so that ergodic- 
ity is restored. In the case of the sandwich, this question 
can be asked in two ways: 

1. How big must the wall separation be so that the 
liquid inside behaves like the bulk? 

2. Given a very large (or infinite) cavity, how far from 
the walls must one look so that the liquid behaves 
like in the bulk? 

The first question implies that one is observing the over- 
lap as a function of d at some reference position within 
the sandwich (typically at the center, since influence of 
the walls very near the interface is always expected). In 
the second question, one considers the overlap as a func- 
tion of z at fixed, very large d. As we shall see, the two 
questions have qualitatively and quantitatively different 
answers. 



A. Point-to-set correlation length 

We first study the decay of the overlap following the 
point-to-set prescription, i.e. measuring the overlap at 
the center of the sandwich (z = d) and varying the dis- 
tance d between the walls (by symmetry, we can actually 
average the overlap over the whole central plane). We call 



FIG. 2. Overlap at the center of the sandwich vs. sand- 
wich half-width d in the parallel setup for (from left to right) 
T = 0.482, 0.350, 0.246, 0.202. Lines are exponential or near 
the centercompressed-exponential fits (see text). Inset: same 
data in semilog plot. 



this point-to-set overlap, computed at the central plane 
q c {d). The behaviour of this quantity is shown in Fig. [2] 
for four different temperatures. The scale of decay of this 
function defines the point-to-set correlation length £. 

A notable feature of q c {d) is that its decay crosses over 
from simple exponential at high temperatures to non- 
exponential at low temperatures. In the low T phase 
a simple exponential fit does a very bad job, while the 
curves can be fitted via a "compressed exponential" form, 

(?c (d) = r!exph(d/0 <: ]+go , (2) 

where the anomaly exponent ( measures the deviation 
from exponcntiality. This specific form is by no means 
the only one capable of capturing the nonexponcntial 
shape. The relevant point is that such nonexponcntial 
behaviour is present, and that it is useful to have a scalar 
parameter (in this case £) to quantify it. 

At high temperatures a semilog plot shows that the 
curves are reasonably exponential, so in order to avoid 
ovcrfitting we fix £ = I and we fit the data to a pure 
exponential. On the other hand, at low temperatures 
there is a clear deviation from exponcntiality (inset of 
Fig. [2]), so that the nonexponcntial fit (Eq.[2|) is used. At 
the lowest T we obtain ( = 2.7 ± 0.2 (see Table I for all 
values of Q. 

This progressive sharpening of the decay at low tem- 
peratures (growing of the anomaly exponent) is also 
found in the spherical cavity [l9j |. but the numerical 
value of the exponent £ is different (lower) in the sand- 
wich case. Thus the geometry of the system may in- 
fluence the strength of the exponential/non-exponential 
crossover but the existence of the crossover itself seems 
not to depend on the geometry and it is therefore a robust 
result. 
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FIG. 3. Overalp in the paralell setup as a function of z for 
a sandwich of half-width d — 4 at several temperatures. The 
free region is wide enough that the overlap can reach its bulk 
value go near the center. 



FIG. 4. Overlap vs. d — z (distance from wall) at fixed d (same 
data as Fig. [3} with pure exponential fits. Inset: same data 
in semilog scale. 



IV. MOSAIC IN THE SANDWICH 



B. Penetration length 

We now consider the decay of the overlap off one single 
wall. It is clear from Fig. [3] that for a large enough values 
of the sandwich width d, the overlap has enough room 
to decay to its liquid value qo at the central plane, at 
all temperatures. Therefore, the decay of the overlap, 
q(z, d ;§> £), as a function of the distance z from one of 
the two walls, is perfectly equivalent to the decay of the 
overlap from a single wall in a semi-infinite geometry. We 
call this quantity simply q(z). 

Fig. [4] shows the behaviour of q(z) focusing on one sin- 
gle wall. The first feature that we notice, in comparison 
with the point-to-set correlation function, is that at all 
temperatures data are well fitted with a simple exponen- 
tial, 



q{z) = exp[-z/X] + q Q , (3) 



where A is the penetration length. Independently from 
the fit quality, the pure exponential behaviour is evident 
from the semilog plot (inset of Fig. @| . This result is in 
agreement with the results obtained with a single wall in 
[20L [2ll | : an exponential decay at all temperatures with no 
sign of crossover to non-exponentiality. This feature is a 
remarkable difference with respect to the nonexponential 
point-to-set correlation q c (d). Such difference is perhaps 
not surprising: the two quantities are conceptually not 
the same, as we shall argue in the next section. 

Apart from the functional form of the decay, another 
difference that we immediately notice is that the pene- 
tration length A seems to be smaller than the correlation 
length £. We will return to this in Sec. I VII 



A. Naive argument 

We have shown in the previous subsection that the 
decay of the overlap in the sandwich has the same 
cxponcntial/super-cxponcntial crossover with tempera- 
ture that is observed in the spherical geometry. Such 
anomalous non-exponential behaviour at low T was ex- 
plained in by using a generalization of the RFOT 
framework. In this Section we will show that, at least at 
the naive level, the same RFOT arguments that hold in 
the cavity can be also applied to the present sandwich 
geometry 

The basic idea of RFOT is that the relaxation of a con- 
fined system is regulated by trade-off between a cost and 
a gain of exploring states other from the one fixed into 
the amorphous boundary conditions. The cost is the free 
energy the system has to pay to form an interface when 
it changes state, whereas the gain is the entropic surplus 
the system enjoys by changing state (28M30| . The slight 
complication of the sandwich is that one must careful to 
take account of its anisotropic geometry. Unless we are at 
some very specific value of the parameters (that we shall 
discuss later) , it seems reasonable to assume that the re- 
arrangement of the mobile part of the sandwich happens 
independently within uncorrelated regions, whose longi- 
tudinal size is larger that the correlation length £. If we 
call A and B two such regions, we are saying that 

Za+b ~ Z A Z B . (4) 

This means that the overlap of the mobile particles will 
be basically a longitudinal average of the overlaps of such 
uncorrelated regions, 

n 

qc (d,L) = -^qi r \d), (5) 
r=l 

where n is the number of uncorrelated regions along the 
sandwich. 
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If we accept this, then the RFOT argument can be run 
over one independent region of longitudinal size ~ £, and 
of width <~ d. Exactly as in the cavity, the entropic gain 
is (all relations are given in the three dimensional case) 

AF gain ~ T^i 2 d. (6) 

The surface tension cost, however, is trickier than in the 
cavity. On one hand, we know that it must scale like 
a length to the power 9, the stiffness exponent. On the 
other hand, we also expect from extensivity reasons that 
this cost must scale like the longitudinal size of the rear- 
ranging region to the power d — 1 = 2: surely, if we build 
a super-sandwich by putting many sandwiches aside, the 
total cost must be additive. We can encapsulate these 
two requirements by writing 

AF cost ~ Yd e f(d/£), (7) 

where f{d/£) is a scaling function that, due to extensiv- 
ity, must obey the relation 

fWO ~ e/d 2 , £ » 1. (8) 
In the end, we get 

AF cost ~ Y?dP- 2 . (9) 

As usual in the RFOT argument, we obtain the corre- 
lation length, i.e. the lengthscale at which the overlap 
decays to zero, as the value of d where the two contribu- 
tions balance, Af ga j„ ~ AF cost . This yields 

dRFOT ~ i^J • (10) 

This is the same prediction as RFOT gives in a cavity 
geometry. This sharp RFOT scenario should then be 
smoothed by including the surface tension fluctuations, 
following Ref. [l9j]. In this way one gets a q c (d) that 
decays on a scale c?rfot, and whose decay is sharper 
and sharper (larger exponent £) the lower the temper- 
ature, in agreement with what we find numerically. In 
this context, the point-to-set correlation length £ must 
be identified with the RFOT lengthscale i^rfot j 

£ ~ dRFOT- (11) 

Note that, at the level of this naive treatment, the differ- 
ence between the two walls vs. the single wall geometry, 
and therefore the difference between q c (d) and q{z), is 
quite clear. In the single wall case the entropic gain is 
infinite, as flipping the entire semi-plane is an advan- 
tage over any interface energy. So, we do not expect 
any trade-off in that case. However, even in the single 
wall geometry, the best distance A where to locate the 
interface will be nontrivial, since it may be entropically 
inconvenient for the system to squeeze the interface too 
close to the wall. But it will be the entropy of the rough 
interface, not that of the bulk, to matter. For this rea- 
son, we do not expect the growth of A to be regulated 
by a classic RFOT trade-off, while do we expect so for 
the point-to-set length £. No surprise, then, that the two 
quantities arc different. 



B. Sharpening in the thermodynamic limit? 

As we argued before, an advantage of the sandwich ge- 
ometry over the cavity is that one can tune the width d 
and the longitudinal size L independently. This means 
that (at least in principle) in the sandwich one can per- 
form the thermodynamic limit L — > oo at fixed d. How- 
ever, because of the statistical factorization hypothesis 
(Eq. S]), the longitudinal size L plays no role at all in the 
RFOT argument. In general, this is not necessarily cor- 
rect. It has been argued in Ref. [3l[ that, depending on 
the specific system's geometry and on the dimensionality, 
the limit L — > oo can actually turn the q c (d) smooth de- 
cay with d, into a bona fide, sharp transition at d = (ipxs j 
even at T > Tfc. We will only sketch the argument here. 

What we have disregarded above is the interaction be- 
tween the different rearranging regions in the mobile part. 
Consider two neighbouring regions, A and B, and ask 
which is the propensity of A to decorrelate from its ini- 
tial state. Clearly, this depends on the frozen boundaries 
enclosing A, but also on the state of the neighbouring 
particles in B [H, [33|. The state of B may favour or 
not the ergodization of A [32|, [33| , and one should take 
into account this interaction. According to the theoreti- 
cal scenario of Ref. [3l| , it turns out that exactly at the 
transition point d = c?rfot, this longitudinal interaction 
can make the sandwich long-range correlated along the 
longitudinal plane. This phenomenon would work in the 
direction of making the transition between IN and OUT 
states sharper and sharper. However, in [3lj it is also 
remarked that such transition is smoothed by the pres- 
ence of the disorder (disorder is generated by the surface 
tension fluctuations along the sandwich) and that this 
has the effect to suppress the transition in a d = 3 sand- 
wich, which is our case. Therefore, one should not expect 
any particular effect when increasing L (the transition 
would not be suppressed in a d = 4 sandwich, nor in a 
d = 3 system with randomly frozen particles, though — 
see Ref. H). 

We report the overlap q c (d,L) for two different sizes, 
L = 16 and L — 25, in Fig. [5] Indeed we do not 
find any evidence of a sharpening of the decay of for 
larger L, which confirms the expectation above. For a 3- 
dimensional sandwich, thus, the naive RFOT argument 
provided at the beginning of this Section is probably good 
enough. 



V. "ANTI-PARALLEL" BOUNDARY 
CONDITIONS 

We now turn to the study of the excess energy pro- 
duced by forcing an interface in the mobile part of the 
sandwich. To do this we use "anti-parallel" boundary 
conditions: we freeze particles on one wall in a configu- 
ration a, and those on the other wall in a different config- 
uration j3 (see Fig. [TJ bottom) . The reason to study this 
geometry is twofold. First, the surface free energy cost is 
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FIG. 5. Overlap at center vs. sandwich half-width d at T 
0.203 and two values of L. 



a key ingredient of the RFOT theory [3fJ, [3J], but little 
is known about it. The very possibility of measuring the 
surface tension between amorphous states is at present 
under debate [35| . and recently arguments against the 
existence of amorphous domain walls (a question which 
is closely related to the previous one) across the bulk of 
the glassy liquid have been proposed (33|. Moreover it 
is not clear if the surface tension in a supercooled liquid 
is a purely entropic phenomenon or if it also includes an 
energy part due to the mismatch of different states. 

Secondly, the numerical study of the excess energy pro- 
vides in principle a method to estimate the stiffness ex- 
ponent another crucial player in the RFOT formulas, 
regulating the growth of the correlation length. Unfortu- 
nately, we shall see that, although the sandwich geometry 
is in principle ideal to determine through the technique 
of the aspect ratio scaling, in practice the present values 
of the correlation length are not large enough to make an 
unambiguous determination of 9. 



A. Interface energy 

We define the excess energy as the difference between 
the extensive energy of the mobile part of the sandwich 
with "anti-parallel" boundary conditions and "parallel" 
boundary conditions, 



AE(d) = E a p(d) — E aa , 



(12) 



averaged over 16 samples. It is important to note that, 
if we average over a large enough number of samples, the 
energy E aa is equal to the (extensive) equilibrium energy. 
In the anti-parallel case, relaxation of the energy is in 
general very slow, and it get slower for smaller d. This 
fact is true even using the accelerated swap dynamics 
(Fig. [5]). Clearly, the system is unhappy with the a/3 
boundary, likely because of the forcing of a domain wall. 
For this reason, at the smallest values of d we do not 
reach a plateau of the energy even for our longest time. 
In these cases we extrapolate the limiting value of the 
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FIG. 6. Excess energy per mobile particle at T = 0.246 and 
several values of d, together with the bulk (PBC) average 
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FIG. 7. Excess energy vs. d at several temperatures and 
L = 16. Lines are exponential fits. Inset: same data in 
semilog scale. 



excess energy by using a power-law fit, 

AE(t,d) = AE(d) + At- a . 



(13) 



Fig. [7] shows the excess energy AE(d) for all the values 
of d studied. As expected, AE decays when increasing 
d, and, at fixed d, the excess energy grows upon lowering 
the temperature. At the largest temperature T = 0.482 
though, AE is basically always zero except for the small- 
est d. The fact that the energy cost to match independent 
amorphous configurations vanishes for high temperatures 
seems to support the existence of the spinodal crossover 
proposed in [36T - l38j . 

At the three lowest temperatures the excess energy 
seems to be well described by an exponential decay with 
d, 



AE(d) 



-d/it 



(14) 



(see inset of Fig. [7]). We must note that, at variance 
with the case of the point-to-set correlation, in the excess 
energy wc do not find any hint of non-cxponcntiality. 
Moreover, the decay of AE(d) defines a new lengthscale 
I. We shall investigate in the next Section whether I can 
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be identified with the point-to-set correlation length £ or 
with the penetration length A. 

B. Aspect ratio scaling 

It is interesting, and potentially useful, to analyze the 
excess energy using some simple scaling relations, par- 
tially inspired by the aspect-ratio-scaling technique in- 
troduced in ref. [39|. The basic ideas of this subsection 
have been already used in the naive RFOT argument of 
the previous section. 

The relevant lengthscales for AE are d, I and L, the 
longitudinal size of the sandwich. The first thing we can 
say is that the excess energy will scale like a length to the 
exponent 9 (which is basically a definition of the stiffness 
exponent). Hence 

AE ~ YL 9 f(d/ L, 1 1 L), (15) 

where Y is the (generalized) surface tension. One can 
choose any of the three lengths to fix the dimensions by 
appropriately changing the scaling function /. The sec- 
ond requirement is that the energy and the excess energy 
must be extensive: in the limit L>ci, the AE from dif- 
ferent pieces of the surface must add up. This implies 
that, 

AE~L 2 , (16) 

a relation very well obeyed by our data. For this to be 
true we need that 

f(d/L 7 d/l)~(L/df- 9 g(d/l). (17) 

Moreover, as we have seen from the data, the scale I 
seems well set by an exponential decay, so it is reasonable 
to assume g(x) = e~ x , so that 

AE(d) ~ YL 2 Jzee~ d/l - (18) 

This is an interesting formula, and one could in principle 
use it to fit the stiffness exponent 9. In particular, the 
formula suggests that, if a purely exponential fit is sat- 
isfactory (as in our case), then 9^2. In practice, the 
formula is useful to discriminate different values of 9 only 
for large d; but for AE to be nonzero at large d, we need 
very large values of /, i.e. very low temperatures, which 
we do not have. In fact, any exponent 9 in the inter- 
val [1, 2] does an equally good job in fitting our data for 
AE(d). In particular, distinguishing between 9 = 3/2 
and 9 = 2 is completely out of the question. Yet, the 
method is conceptually interesting, and future simula- 
tions, at lower T, may eventually use it to determine the 
stiffness exponent. 

VI. COMPARISON OF THE DIFFERENT 
LENGTHSCALES AND THE ISING CASE 

Let us summarize the three lengthscales we have mea- 
sured. The first one is the point-to-set correlation length, 



£, defined as the decay scale of the overlap measured at 
the centre of a sandwich of half-width d, 

q c (d)~exp[-(d/0 C ]+<Zo. (19) 

Previous investigations suggest that this is the true static 
correlation length of the system, the one relevant for the 
structural rearrangement [40j . Moreover, there is evi- 
dence [H that £ has to be identified with the RFOT 
correlation length, discussed above. The remarkable fea- 
ture of the point-to-set correlation length is that its asso- 
ciated correlation function has a non-exponential decay 
at low temperature. 

The second lengthscale is the penetration length A, 
which regulates the decay of the overlap off a single wall 
in a semi-infinite geometry, 

q(z) ~cxp[-z/A] +q . (20) 

This lengthscale seems to have a different physical mean- 
ing than the correlation length £, as also suggested in [3l| . 
It seems to embody the extent to which the effect of a 
single frozen wall penetrates into the system, rather than 
the average size of a rearranging region. At variance with 
the point-to-set correlation length, A regulates a purely 
exponential decay of the overlap. 

Finally, we measured the lengthscale I associated to 
the decay in d of the excess energy produced by imposing 
"anti-parallel" boundary conditions, 

AE(d) = E a0 - E aa ~ exp[-d//]. (21) 

As in the case of the penetration length, the excess energy 
lengthscale I is associated to a purely exponential decay, 
at least down to our lowest available temperature. 

What can be said about the quantitative relationship 
(if any) between these three lengthscales? We report 
them all in Table IVI| together with the anomaly expo- 
nent £. 

TABLE I. Point-to-set correlation length £, and anomaly ex- 
ponent £ , from a fit of Eq. 1191 penetration length A, from a 
fit of Eg. 1201 and excess energy decay lengthscale I, from a fit 
of Eq. [21] At the highest temperature the value of I has large 
uncertainty as we have very few nonzero values of AE. 
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One could object that much of the comparison depends 
on the fitting procedure of the data, which is not a nice 
thing. This is certainly a concern. However, we notice 
that extracting the lengthscales by crossing the various 
functions with arbitrary threshold would not be any bet- 
ter, for two reasons: first, in presence of a nonexponcntial 



2 - 
1 . 8- 
1 . 6- 
1 . 4- 
1.2- 

1 - 
0.8- 
. 6- 
0.4- 



0.4 0.6 0.8 



1.2 1.4 1.6 1.1 



10 
8 

6 - 
4 
2 



10 



FIG. 8. Penetration length A vs. point-to-set correlation FIG. 10. Energy decay length vs. two-spin correlation length 

length f f° r the Ising model in the square lattice. Data are from 

Monte Carlo simulations on a 100 x 100 lattice with single-flip 
Metropolis dynamics performed above the critical point, at 

2 1 temperatures T = 2.5 J, T = 2.4J, T = 2.35J and T = 2.32 J, 

1 • ^~ where J is the Ising coupling constant (the critical point is 

1 - 6- T c m 2. 269 J). The correlation length was obtained from a fit 

1 . 4- of the spin-spin space correlation function C(r) = (S(0)S(r)). 

1 . 2- To determine the length I, sandwich configurations were pre- 

1-1 i - pared as explained for the liquid case, measuring the excess 

. 8 - ..---'*' energy AE(d) = E a p - E aa for d = 1, 2, 3, 4, 5, 7.5, 10, 15, 

q g 20, 25 and 45 lattice spacings, and fitting to an exponential 

; 4 : m . ■ decay. 
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FIG. 9. Energy decay length I vs. point-to-set correlation 
length f . 



decay (as q c (d) unmistakably has), with a T-dcpcndcnt 
anomalous exponent £, the arbitrary value of the thresh- 
old can strongly bias the dependence of £ on T; secondly, 
these are dimensionally different, inhomogeneous func- 
tions, so it would be hard to choose coherently a crossing 
point for each of them. An honest fit is the best we can 
do. 

From the table we see that the correlation length £ is 
larger than the other two scales. Of course, what really 
matters is their mutual T-dependence, namely: is there 
any of them that grows significantly faster than the oth- 
ers, when lowering T? More precisely, we would like to 
understand whether or not these lengths are ruled by dif- 
ferent exponents. Because of this, comparing the three 
plots, £(T),\(T),l(T), is not a good idea: constant fac- 
tors would show up as increasing differences, conveying 
the (wrong) idea that one length is growing faster than 
the other. The best thing to do is to plot one lengthscale 
vs. the other, parametrically in T. This is what we do in 
Figs. U and [S 

The result of this comparison is unfortunately not con- 
clusive. Even though, as we already said, the correlation 
length £ is larger than the other two, all mutual depen- 



dences are not far from linear. This means that, with 
such data, we cannot claim that £ is growing with an ex- 
ponent significantly different from the other two, which 
would be the only proof of a qualitative difference be- 
tween these scales. Of course, our data do not either rule 
this out. 

In such a murky situation, some conceptual help may 
perhaps come from the well-known Ising model. There 
one can use (true) anti-parallel boundary conditions to 
force a domain wall (below the critical temperature Tc). 
Then there is a finite surface tension and the excess free 
energy (anti-parallel minus parallel) in the limit d — > oo 
tends to the finite value aL d ~ 1 1 where a is the surface 
tension. For T > T c , one can instead impose (similarly 
to what we have done above) two different paramagnetic 
configurations on the two sides of the sandwich and mea- 
sure the excess energy. In this case we expect the excess 
energy to decay to zero for large d, but on what scale does 
this happen? Simulations in two-dimensions show that 
the excess energy decays exponentially with a lengthscale 
I , which we can compare with the Ising correlation length 
£, calculated from the standard spin-spin space correla- 
tion function. Data show that I and £ scale linearly with 
each other (parametrically in T, Fig. [TO]) , though I is 
somewhat 2 smaller than £. Therefore, in Ising, correla- 
tion decay and excess energy decay seem to track each 
other quite closely. 
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VII. CONCLUSIONS 

Similarly to what we found in the spherical geometry 
[l9j . the sandwich data show a crossover from exponen- 
tial to nonexponential relaxation of the point-to-set cor- 
relation function, upon lowering the temperature. We 
remark that such crossover is one of the very few (if not 
the only one) static landmarks differentiating at the qual- 
itative level the fluid phase from the deeply supercooled 
phase in glass-forming liquids. Having found this feature 
now in two different geometries makes it quite a robust 
phenomenon. 

Up to now, the only reasonable explanation of such 
sharper-than-exponential relaxation has been given in 
the context of RFOT. We reported a naive RFOT ar- 
gument for the sandwich and showed that there should 
be no essential variations (at least in 3 dimensions) with 
respect to the standard argument one uses in the cavity 
geometry. Our sandwich results therefore give further 
support to the theoretical connection between point-to- 
set nonexponential relaxation and RFOT. 

A more thorough formulation of the RFOT argument, 
based on a renormalization group framework, could turn 
the nonexponential, but smooth, drop of the overlap at 
d = g^rfot j into a true transition, in the limit L — > oo 
[3l| . The very existence of such limit would be one of 
the main benefits of the sandwich vs. the cavity geome- 
try. However, this transition is supposed to be smeared 
out by disorder (surface tension fluctuations) in a three- 
dimensional sandwich j3lj . and indeed, by substantially 
increasing L, we do not find any relevant change in the 
point-to-set correlation function. 

A different experiment consists in measuring the over- 
lap decay off a single all. In this case, we found a be- 
haviour rather different from the point-to-sct correla- 
tion function. First, and most important, this decay is 
purely exponential, even a the lowest temperature stud- 
ied, where the point-to-set correlation function is clearly 
nonexponential. Hence, the single wall seems less than 
ideal to characterize at the qualitative level the deeply 
supercooled phase. 

Second, the lengthscale of this single-wall decay, i.e. 
the penetration length A, seems to be smaller than the 
point-to-sct correlation length, £. According to ref. (4l| . 
both £ and A should diverge at Tk, but with different 
exponents, in particular the divergence of £ should be 
sharper than that of A. This is due to the fact that A is 
not controlled by the RFOT entropy vs. surface tension 
competition mechanism so directly as £ is. Even though 
we are far (to say the least) from the T ~ Tk region 
where the RG arguments of ref. [4l[ hold, we can at least 



say that our numerical data are not in contradiction with 
this scenario. 

We used "anti-parallel" boundary conditions in the 
sandwich to measure the excess energy associated to an 
interface. This quantity seems to decay purely expo- 
nentially with the half-width d of the sandwich, over a 
lengthscale I that grows by lowering the temperature. 
Unfortunately, the value of I we obtain even at the lowest 
T is not large enough to make it possible an estimate of 
the stiffness exponent 8 using aspect ratio scaling. How- 
ever, the technique seems promising in this context, and 
perhaps future simulation will reach a regime able to dis- 
criminate between different (theoretical) values of 6. 

An obvious question is whether and how the length- 
scale of the excess energy / is related to the other two 
lengthscales, and in particular to the correlation length, 
£. We do not have a final answer to this question. Even 
though £ seems to be quantitatively larger than I (about 
a factor 4), there is no clear evidence of a nonlinear con- 
nection between the two lengths. On one hand, by fol- 
lowing an economy criterion, we are tempted to conclude 
that point-to-set and excess energy are regulated by one 
lengthscale, as it happens in the Ising model above T c . 
On the other hand, the very different kind of relaxation 
(nonexponential for the point-to-set correlation function, 
exponential for the excess energy), and the fact that I is 
significantly smaller than £, seem to suggest otherwise. 
We cannot but leave the question open. 

Finally, the Ising example calls for caution in the inter- 
pretation of our results about the excess energy. The fact 
that in a purely paramagnetic state one finds a behaviour 
of AE vs. d so similar to the glass-forming case, means 
that a finite AE is not by itself proof of the existence of 
a surface tension. This calls for a through investigation 
of the possible entropic contribution to the glassy surface 
tension (which we do not measure here). 
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